Neutron spectroscopy and magnetic relaxation of the Mn 6 nanomagnets. 
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Inelastic neutron scattering has been used to determine the microscopic Hamiltonian describing 
two high-spin variants of the high-anisotropy Mn6 nanomagnet. The energy spectrum of both 
systems is characterized by the presence of several excited total-spin multiplets partially overlapping 
the S=12 ground multiplet. This implies that the relaxation processes of these molecules are different 
from those occurring in prototype giant-spin nanomagnets. In particular, we show that both the 
height of the energy barrier and resonant tunnelling processes are greatly influenced by low-lying 
excited total-spin multiplets. 

PACS numbers: 75.50.Xx, 75.40.Gb, 75.60.Jk, 78.70.Nx 
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I. INTRODUCTION 

Molecular nanomagnets (MNMs) have recently at- 
tracted considerable interest because at low temperature 
T they display slow relaxation of the magnetization M of 
purely molecular origin [l[ . The main relaxation mecha- 
nism is provided by the interactions of the spin degrees 
of freedom with phonons, either through modulation of 
the local crystal fields on individual magnetic ions, or 
through modulation of two-ion interactions. Typically, 
the relaxation dynamics are modelled by restricting the 
spin Hilbert space to the ground total-spin multiplet 
only, as resulting from the usually dominating isotropic 
exchange interactions. In this framework, each ./V-spins 
molecule is described as a single giant spin S in an 
effective crystal-field potential. At temperatures of 
few K, the reversal of M occurs through a multi-step 
Orbach process yielding a thermally activated behavior 
of the relaxation time, r = T$exp(U / ksT), where the 
energy barrier U is set by the effective axial anisotropy 
experienced by the giant spinQ. 

We show that excited S'-multiplets strongly influence 
the energy barrier for the relaxation of M in two 
high-spin (S = 12) variants of the high-anisotropy Mn6 
nanomagnet (a record barrier U — 86.4 K, and U = 53.1 
K [!, 0]). I n fact, inelastic neutron scattering (INS) and 
Frequency Domain Magnetic Resonance Spectroscopy 



(FDMRS) show that the two variants are characterized 
by a similar anisotropy but have significantly different 
exchange interactions. The large difference in U is 
mainly due to the presence of relaxation paths passing 
through excited S multiplets partially nested within 
the ground one. In addition, because of ^-mixing in 
the wavefunctions, these excited manifolds may lead to 
resonant inter-multiplet tunnelling processes for fields 
of a few thousands of Gauss. These are associated with 
additional steps in hysteresis cycles which are absent in 
the giant spin model. 



II. SPIN HAMILTONIAN AND MOLECULAR 
ENERGY LEVELS 

The two Mn6 molecules have chemical formula 
Mn 6 2 (Et-sao) 6 (0 2 CPh(Me)2)2(EtOH) 6 (higher bar- 
rier) and Mn 6 2 (Et-sao) 6 (0 2 CPh) 2 (EtOH) 4 (H 2 0)2 
(lower barrier) and are nearly isostructural. The mag- 
netic core comprises six Mn 3+ ions arranged on two 
triangles bridged by oxygen atoms (Fig. 1). Each Mn 3+ 
ion has a distorted octahedral cage of ligands, with the 
Jahn— Teller axes all approximately perpendicular to the 
planes of the triangles. The six Mn 3+ ions have spin s 
= 2 and are coupled by dominant ferromagnetic interac- 
tions, leading to a high S = 12 total-spin ground state, 
as can be inferred by magnetization measurements, see 
Fig. 2 Hi- 

Each Mn6 molecule can be described by the following 
spin Hamiltonian Q 

H = J ijS (i) • sfj) + disl{i)+ 

i<j i 
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FIG. 1: (color online) Structure of the high-barrier Mn6 
molecule (H omitted for clarity). Green : Mn, red : O, blue : 
N, dark grey : C). 
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FIG. 2: (color online) Measured (points) and calculated 
(lines) T-dependence of the susceptibility for the two Mne 
molecules. Calculations are made by including exchange in- 
teractions only, with the parameters determined by INS. g = 2 
has been assumed. 
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where s(i) are spin operators of the i th Mn ion and 
S = J^i s (*) is the molecule's total spin. The first term is 
the isotropic exchange, while the second and third terms 
describe axial local crystal-fields (a z axis perpendicular 
to the triangles plane is assumed), and B is the external 
field H' (neglected in the following) represents addi- 
tional small anisotropic terms. The minimal set of free 
parameters is given by three different exchange constants 

JlV = J%, Jl2 = J23 = J 13 — Jl'T — J2'3' = Jl'3' = J2, 

and Ji 3 ' = Ji'3 = J 3 (Fig.l) and two sets of crystal- 
field (CF) parameters d\ — dy, c\ = ci<, and di = d^ 1 
c 2 = c 2 ' . The ligand cages of sites 1 and 3 are rather sim- 
ilar and we assumed the corresponding CF parameters to 
be equal. Since experimental information is insufficient 
to fix independently the two small c parameters, we have 
chosen to constrain the ratio c\jci to the ratio d\jdi. 
The anisotropic terms break rotational invariance and 



here lead to a large amount of mixing of different S mul- 
tiplets (S mixing [6] ) . In the following we label the states 
by their leading S'-component. 

To determine the parameters appearing in (1), we 
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FIG. 3: (color online) INS spectra for the higher barrier 
molecule collected on IN5 with incident wavelength of 3.4 A 
for T= 2 K (black) and T = 16 K (red). Lines are theoretical 
calculations with model (1). 
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FIG. 4: (color online) INS spectra for the higher barrier 
molecule collected on IN5 with incident wavelength of 6.7 A 
for T= 2 K (black) and T = 16 K (red). Lines are theoretical 
calculations with model (1). 

have used the time-of-flight neutron spectrometers NEAT 
at the Hahn Meitner Institut (Berlin) and IN5 at the 
Institute Laue Langevin (Grenoble). In addition, we 
have performed FDMRS measurements at the Univer- 
sity of Stuttgart. Since FDMRS is sensitive to intra- 
multiplct transitions only, its use in conjunction with 
INS makes easier to assess the character of the differ- 
ent observed excitations. Figs. 3 and 4 show exam- 
ples of INS spectra together with theoretical simulations. 



3 



The analysis of INS and FDMRS data leads to the fol- 
lowing parameters for the higher (lower) barrier com- 
pounds : Ji = -0.84(-0.61) meV, J 2 = -0.59(-0.31) 
meV, J 3 = 0.01(0.07) meV, di = -0.20(-0.23) meV, 
d 2 = -0.76(-0.97) meV and c x = -0.001(-0.0008) meV 
[5J. These parameters are consistent with susceptibil- 
ity measurements (see Fig. 2). Hence, anisotropy is 
similar in the two variants whereas the dominant ferro- 
magnetic exchange is substantially larger in the higher- 
barrier compound. Figures 5 and 6 show the energies 
of the S'-multiplets resulting from the exchange part of 
(1). It is evident from these figures that the exchange 
splitting between the ground S = 12 manifold and many 
excited multiplets, including low-S ones, is smaller than 
the energy scale of anisotropic terms in (1). In particu- 
lar, Fig. 6 shows that in the lower-barrier compound the 
energy of the lowest-lying S = multiplet is only about 
4 meV larger than that of the ground one. Thus, the 
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FIG. 6: Energy of all S'-multiplets for the lower barrier com- 
pound, resulting from the exchange part of (1) with the pa- 
rameters given in the text. 
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FIG. 5: Energy of all S'-multiplets for the higher barrier com- 
pound, resulting from the exchange part of (1) with the pa- 
rameters given in the text. 

giant spin mapping completely breaks down in these two 
molecules, not only for the large 5-mixing in the wave- 
functions, but even for failing to account for the number 
of states located below the barrier. The presence of a 
center of inversion (characterizing the structure deter- 
mined at 150 K) implies that exchange multiplets can 
be classified according to their parity with respect to the 
associated spin-permutation operation. For instance, the 
ground S — 12 states are even, whereas the lowest- lying 
5 = 11 states are odd. 



adopt the well established framework presented in [7J, |8[ 
for the irreversible evolution of the density matrix pit). 
By focusing on time scales much longer than the typi- 
cal periods of free evolution of the system, the so-called 
secular approximation enables the time evolution of the 
diagonal matrix elements of the density matrix to be de- 
coupled from that of the off-diagonal ones. In particular, 
the populations of the molecular eigenstates, pi{t), evolve 
through master equations: 

Piit) = ^2w lmPm it), (2) 

m 

where Wi m is the Im element of the rate matrix, i.e., the 
probability per unit time that a transition between levels 
|m) and \l) is induced by the interaction with the phonon 
heat bath. The latter can be calculated by perturbation 
theory once magnetolelastic (ME) interactions have been 
modelled. Experimental information is totally insuffi- 
cient to fix all the many possible parameters appearing in 
the ME coupling potential. By assuming that the main 
contribution to this coupling arises from the modulation 
of the local rank-2 crystal fields and that quadrupolc mo- 
ments of each individual Mn ion are isotropically coupled 
to Debye acoustic phonons, we obtain [9j 

Wi m = 7 2 Af m n(A /m ) ( l \°Qi,iM)Hx 



III. RELAXATION OF THE MAGNETIZATION 

The main difference between the two Mii6 variants is 
the position of the excited S manifolds, therefore the 
comparison between the relaxation behavior of these sys- 
tems offers the opportunity to study the role played by 
low-lying excited multiplets. To address this issue, we 



(m\O quq2 ( Sj )\l), (3) 

where i and j run over Mn ions, 9ig2 (si) are the 
components of the cartesian quadrupole tensor operator, 
n(x) = (e Hx / kBT -l)- 1 and A /m = (Ei-E m )/%. The sin- 
gle free parameter 7 is proportional to the ME coupling 
strength. In spite of the complexity of their energy spec- 
trum, for both variants the resulting relaxation spectrum 
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at low T is dominated by a single relaxation time display- 
ing a nearly Arrhenius behavior r(T) = r ey^p{U / KgT), 
in agreement with experiments [5(. In particular, the 
effective energy barrier for the magnetization reversal is 
crucially dependent on the position of low- lying excited S 
multiplets. Indeed, the large difference between the bar- 
riers of the two molecules mainly results from the vari- 
ation of the ferromagnetic exchange constants and not 
from a variation in anisotropy. It is also worth to note 
that in the present case U is not set by the energy of the 
lowest-lying Ms = state as one could naively expect. 
The existence of several total-spin multiplets partially 
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FIG. 7: Energy levels of the higher-barrier compound as a 
function of B, applied along the z axis. The color maps for 
each state the value of 5 e //, where < S 2 >:= S e ff(S e ff + 1). 
Ellipses and the square indicate representative examples of 
level crossings. 

overlapping with the ground one implies that a magnetic 
field applied along the molecule's easy axis induces cross- 
ings involving states which are absent in a giant spin de- 
scription. These can be turned into anticrossings (ACs) 
by the small transverse terms contained in H' . Figure 7 
shows that besides the "traditional" intra-,? = 12 cross- 
ings at about 0, 0.4 (e.g., the square) and 0.8 T, many 
more occur for intermediate values of B. In particular, 
there are three different kinds of " non-traditional" field- 
induced crossings, i.e., which are absent in the giant spin 
description of the molecule. First of all there are cross- 
ings involving a pair of states belonging to a single low- 
energy S multiplet different from the ground one (for in- 
stance the higher-lying ellipsis at 0.3 T in Fig. 7). These 
crossing may lead to relaxation shortcuts if the domi- 
nant relaxation path passes through them. The other 
two kinds of crossings (e.g. the other ellipses in Fig. 7) 
involve pairs of states belonging to different S'-manifolds. 
The distinction arises from the symmetry properties of 
the two crossing states, i.e., they may have the same par- 



ity or not. In the first case (e.g., the crossings between 
red and green curves in the two ellipses at about 0.6 and 1 
T) H' usually leads to an AC, whereas in the second case 
(e.g., the ellipse at about 0.3 T) an AC may occur to the 
extent that at low T the inversion center is removed by a 
distortion, leading to terms with low-enough symmetry 
in H' . Even if there are no structural data below 150 
K, it is not unlikely for magnetic molecules to undergo a 
symmetry lowering at low T (see, e.g., [l(J EH)- "Non- 
traditional" ACs with the associated resonant incoherent 
tunneling, may result in relaxation shortcuts leading to 
additional steps in hysteresis cycles, absent in a giant- 
spin picture. For instance, Fig. 8 shows an example of 
derivative of the hysteresis curves measured at T = 3 K 
for the higher-barrier molecule Q. The "traditional" ACS 
produce minor features in these curves, and the two main 
peaks are associated with " non-traditional" ACs. The ef- 
fect of low-field ACs is more evident at higher t[||. For 
instance, there are features at 0.3 T which may a priori 
originate from both crossings marked on Fig. 7. 
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FIG. 8: Example of derivative of the hysteresis curves mea- 
sured for the higher-barrier molecule in [3|. For each value 
of B, there are two points corresponding to increasing or de- 
creasing B in the hysteresis cycle. 



IV. CONCLUSIONS 

By exploiting INS and FDMRS we have determined 
the microscopic Hamiltonian of two different high-spin 
variants of the Mn6 nanomagnet. We have found that 
excited S multiplets overlapping with the ground one 
strongly affect the magnetic relaxation process. More- 
over, we have demonstrated the existence of tunnelling 
pathways involving pairs of states belonging to different 
total spin manifolds. Hence, the energy separation be- 
tween the ground and excited multiplets may be a key 
ingredient in determining the relaxation and tunnelling 
dynamics of molecular nanomagnets. 
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